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O Abstract 

^ An extended analytical model for particle dynamics in fields of a highly-nonlinear plasma wake 

^ field (the bubble or blow out regime) is derived. A recently proposed piecewise model (Kostyukov 

et al., New J. Phys., 12 , 045009 (2010)) is generahzed to include a time dependent bubble radius 



OO 



r"| and full field solution in the acceleration direction. Incorporation of the cavity dynamics in the 

Oh 

' model is required to simulate the particle trapping properly. On the other hand, it is shown that 

«... 

c/^ the previously reported piecewise model does not reproduce the formation of a mono energetic 

Q^ peak in the particle spectrum. The mono energetic electron beams are recovered only when the 

Q full longitudinal field gradient is included in the model. 
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When a plasma wake field is driven to nonlinear amplitudes, it takes the form of an 
electronic cavity that is called "bubble" in the case of a short pulse laser driver [Ij or 
"blowout regime" in the case of dense highly energetic particle bunches [21 13]. For the laser 
driven case, the pulse is shorter than the plasma wavelength and fits perfectly into the first 
half of the plasma period. The laser intensity is thought to be high enough that the created 
wake field breaks after its first oscillation. In this regime, the wake field takes the form 
of a distorted spherical cavity from which all electrons are banished and that moves with 
nearly the speed of light through the plasma. In the following we refer this wake-field as 
"the bubble". 

Inside the bubble, an accompanying electron bunch - the so called beam load - is ac- 
celerated until the spherical shape of the bubble breaks down or the electrons get out of 
phase with the wake. The break down of the sphere happens as soon as the driver energy 
is exhausted [4l[5]. The fields inside the bubble can be calculated from Maxwell's equations 
for an empty sphere [6j and a sphere that takes into account the field enhancement at the 
bubble back due to electron sheet crossing [7, 8j. An even simpler model for the fields inside 
the bubble has been proposed recently: a static piecewise model (PWM) [9j. This simpli- 
fied model neglects any spatial or temporal field dependence, but due to its simplicity it 
is possible to calculate amplitude envelope functions to the transverse particle coordinates 
and momenta, and to analyze some self-injection physics. In contrast to recent experiments, 
the PWM misses the production of mono energetic electron beams. Thus, we are forced to 
modify the PWM to a more general case that remains analytically treatable and improves 
the self-injection physics found in the PWM. 

Our first approach is to consider fields that are time dependent, but piecewise constant 
in space. In this case our model is a sphere with an adiabatically growing radius R{t) and 
a growing trapping cross-section. A physical justification for this approach is that now the 
interaction of the electron bunch with the bubble border is taken into account. Since the 
fields remain constant we still miss the production of mono energetic electron beams. Thus, 
in a second step we introduce a model that treats the right field gradient in the longitudinal 
direction. For this case indeed we find energy spectra that are comparable to PIC simulations 
and a full gradient model. 

Our work is based on the PWM. Thus, similar to [9], our bubble models neglect the 
interaction between the accelerated electrons and the laser pulse at the bubble front. We 



also assume a continuous plasma background density, and a reduced electron dynamics in 
the z-p-plane, where p is the perpendicular distance of an electron to the z-axis and z the 
propagation direction of the laser pulse. The motivation for all these assumptions are several 
observations in the phenomenological model of the bubble regime [6l [TO]. 

The full gradient potentials inside a bubble are given in jHl [TT] . In normalized variables 
they are A\\ = —(f = $/2, and A^ = 0, where $(r) = |rp/4 is the bubble potential. 
The normalization is done by r ^ rcj^g/c for length, t -^ tujpe for time, v -^ v/c for 
velocity, p -^ p/(mec) for the kinetic momentum, and $ -^ e$/(meC^) for potentials. Here 
(jjpe = (47re^no/me) is the electron plasma frequency, rrie is the electron mass, c is the 
velocity of light, and no is the plasma background density. The equations of motion of 
the full gradient model in a co-moving frame of reference are Eq.(l)-(2) in Ref.[10j. In the 
PWM, ^ = z — Vt and p are replaced by 

C = sgn(C)F, p = sgn(p)F, (1) 

where F is the field-strength parameter of the PWM and V the group velocity of the laser 
pulse. 

The time-dependent PWM (tPWM) 

One possibility to examine the interaction between the electron bunch and the bubble border 
is to consider an adiabatically slow bubble growth. This growth significantly affects the 
electron self-injection physics because the bubble back moves slower than the bubble [HI [HI 

[iniiiMS]. 

To adopt the idea of a growing bubble to the tPWM we make the same consideration as 
in [131 E] • Applied to the full gradient model [6j , we find the sufficient trapping condition 
lim^^oo ^{t) = — oc. Thus, a new approach for our wake field potential in the tPWM is 

*(r,t) = f(ri + H)-^, (2) 



with R{t) = (l + 6t)i?o, C = z — V t^V = V — eRo^ V = \J\ — 7o ^, and the gamma-factor 
of the bubble front 70. In these coordinates the bubble center rests and the equations of 
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FIG. 1. (la): Trajectory in the PWM, particles is not trapped, (lb): Trajectory in the tPWM, 
particle is trapped. Both simulations are done with i?o = 5, F = 4, po = 4.92, and 70 = 10. The 
black solid circle has radius Rq. The blue dotted circle is the bubble after the simulation ends. 



(Ic): Simulation configuration: Bubble running into a grid of electrons. 



motion for the tPWM are 
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with 0- = (1 + V-) sgn{^-)F + (1 + et)Rle. 
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Based on the tPWM we perform two kinds of simulation. In the first part we solve the 
equations of motion ([3)l5|) for a single electron and discuss its trajectory, while we always 



assume the initial conditions {^qY + pi = i?o, Po = 0, and e,^o^,po > 0. In the second 
part we place a grid of particles with zero momentum in front of the bubble as shown in 
Fig{Tc]and solve the equations of motion for every individual particle. The resulting energy 
spectra of trapped particles then give us a better insight in the self injection dynamics and 
acceleration mechanism in the tPWM. 



If we compare a particle's trajectory in both the PWM (Fig, la) and the tPWM (Fig, lb), 
we discover a new kind of trapping dynamic. The particle in the static model is lost because 
its trajectory crossed the bubble border. In the time-dependent model, this situation can 
be prevented as the border moves out. The amplitude function of the oscillation in the 






t t iogiBo) 

(a) (b) (c) 

FIG. 2. Envelope approximation functions (red dotted lines) and numerical solutions (blue solid 
lines) to p(t) and p±{t) in the tPWM. The simulation parameters are F = 4, i?o = 6, 70 = 10, 
and e = 0.0004. ( [2c| ): Trapping cross-sections from analytical prediction (blue line) and numerical 
simulations (red circles) for 70 = 10, e = 0.0002, and F = 6 in the tPWM. 



trajectory determines whether the electron will be trapped. So, Kostyukov et al. [9j describe 
a way to calculate the envelope approximation functions p{t) and p±{t) to the amplitude 
functions of p(t) and p±{t). We follow their steps in the limit eR^ <C 2F and introduce the 
times tn by p{tn) = for all n > 1. Then, with n ^ 1 and a guiding center approximation 
to P||(t), we obtain the equation system 



P\\{t) = -[{1 + V-)F-Rle]t, 



~ 2/3 



p^{t) = 1.00c?/o 



p{t) = 4.01- 



_p5/3 



\n-it)\ 



l + V-)F-Rle]'/H'/', 



-1/3 



(6) 
(7) 

(8) 



T2-(t)2[(l + y-)F-i?2e]i/3 

while we use the constant c = c{F, Rq, 70, e, ^q) to fit the curves to our simulations. 

Two exemplary plots of the envelope functions p± (t) and p± (t) (red dotted lines) together 
with the simulated coordinate and momenta of a particle (blue solid lines) are shown in FigJ2J 
The analytical predictions are in good agreement with the simulations and converge very 
close to the actual enveloping functions. 

To calculate the trapping conditions, we demand |rp = (^~)^ + p^ < R? and show that 

for n ^ 1 

Ro^lF 



io< 



+ 



i?nF3/4c3/2 



F ' {F - RlefH" 



1.48 + 0.72 



2F - Rle 



(9) 



is the condition that must hold to trap electrons. From Eq.O) and p^ = R{t) — {^^ ) we 
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FIG. 3. (3a): Energy distribution for F = 4, i?o = 4, e = 0.004, and 70 = 5 in the tPWM. (3b): 



PIC energy spectra taken from Fig. 2 in ^16j. (3c): Energy spectra in the hPWM for 70 = 15, 
i?o = 8, and F = 4. 



calculate the trapping cross section a = ti{R^ — p^j^) with pmin = minpo- If we assume 
Ai? = R — pmin = const^ we find that a{t) is linearly time-dependent. To confirm our 



estimation, we plot log{a) and log(i?o) into one diagram (see Fig, 2c). As we observe, the 
analytical prediction (blue line) follows the numerically found trapping cross sections (red 
circles) and stays especially for large Rq in a 5% error range. 



To analyze the energy spectra in the tPWM a grid of electrons with zero momentum is 



placed in front of the bubble (see Figlc). Afterward, we solve the equations of motion for 
every individual particle, i.e. we neglect the interaction of the test electrons. The electron 
spectra we receive after the electron bunch of trapped electrons has passed its acceleration 



length look like that in Fig, 3a Here, we don't observe any mono energetic electron beams 
but a quasi equal distribution of the particles' energy. Apparently, this contradicts pic 



simulations results as shown in Fig, 3b and experiments [f7l-[T9]. 

The monoenergetic spectra can be recovered in the full gradient model. Thus, in the 
following section we introduce another extension of the PWM that is still analytically treat- 
able and preserves a production of mono energetic electron beams. 



The half-gradient PWM (hPWM) 

In the hPWM we still replace the perpendicular field E± by the constant field-strength 
parameter F. In longitudinal direction, however, we take the full field gradient. The resulting 
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wake-field potential is 

Here, the coordinate ^ = z — Vt always contains the group velocity of the generating laser 
pulse at the bubble front because the bubble radius is fixed. 

If we repeat the simulations of the electron grid in the the hPWM, we get energy spectra 
similar to those shown in FigJ3| As we see, both PIC-simulated spectra in Figj3b] and our 



spectrum from a simulation with F = 4, i?o = 8, and 70 = 15 exhibit the same structure. 
Since the hPWM is able to describe the production of monoenergetic electron beams, we 
calculate the equations of motion 

dP± _ fi^P}[] ^g^(/^)-^ (12) 



dt V 7/ 4 

dt~ J ^' rft " 7 ' ^ ' 

with ri = (1 + V)^/4:. A simple form for the envelope approximation functions to p{t) and 
p±{t) can be found, if we express them in terms of r = t — t^ax^ where ^{tmax) = 0. For 
the calculation we follow [9j again and assume 7(t) <C 2(ncr/^o)^ for all t for the electron 
energy. Then we have 

P\\{r) ^ Pmax - \^ ^ T^ (14) 

^^(r) = VA(Fp||(r))V3, (I5) 

p(r)=A(Fp^||(r))-V3, (16) 

with A = '2}l'^C{tj^^ I^Vlitj^^ . As in the previous section, the constant C(t^) is used to fit 
our approximations to numerical simulations. 

To compare the approximated functions to the actual amplitude envelope functions, we 
solve the equations of motion for a single electron with initial conditions Co + Po = ^o^ 



Po = 0, and <^05Po > 0. An exemplary plot of Eqs.(15 16) together with numerical solutions 
is shown in Figjij The analytical predictions (red dotted lines) are in good agreement to the 
simulations (black solid lines) and converge very close to the actual enveloping functions. 
Similar to Eqs.Q the self-injection physics is described by the limiting condition 

p^i„-F+|4 + F -_--^ + ^_-^l (17) 
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FIG. 4. envelope approximation functions (red dotted lines) and numerical solutions (blue solid 
lines) to p{t) and p±{t) in the hPWM. The simulation parameters are F = b, Rq = 10, and 70 = 20. 



together with 



2AF2/3 1/3 7^2^^-2/3 

P.' + — — z = 0. 



p5/3 _ 7gAi^-/- pi/3 , % 



3 " 3 

?2 ^2 



(18) 
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The trapping cross-section a = 7t{RI — p^i^) i^ constant in time since Rq and pmin are. If we 



plot a and Rq into one diagram (see Fig, 5a), we observe that for a large range of parameters 
the accuracy remains in a 1% range. 

To get some analytical expressions for the limits of the bubble parameters, the so called 



trapping conditions, we start again with Eq.(17) and demand pmin < ^o- Then we find after 
some algebraic transformations 



70^ 



< 2Po(2Fi?o - i?o + 4 - 2AF2/^Po~^^^ + A^P"^/^Po~^^^)"^ 



as well as 



p5/3 



Kq 



R.P}" 



Rq Pq 



p2/3 



A^ 

d2/3 



<0. 



(19) 



(20) 



These two equations give us all parameter limits for the hPWM. To demonstrate that these 
complicated-looking limits are indeed well calculated we plot the largest F for which still an 
electron is trapped against Rq in a simulation with one electron (compare former sections) 
in FigjSb} As we see the analytical found F^ax (blue line) fits perfect to our simulations 
(red dots). 

Conclusion 
Summarizing our expansions of the original PWM we note two things. First, we see that 
it is possible to treat an enhanced interaction between the electron bunch and the bubble 
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FIG. 5. (5a): Trapping cross-section for F — Ro/2. ( |5b| ): Maximal applicable F. Both figures are 
simulated with 70 = 15 and belong to the hPWM. 



potentials. This in turn leads to a growing trapping cross-section and thus to a more efficient 
self-injection. 

Energy spectra in our simulations show that no mono energetic electron beams are gener- 
ated, if all field gradients are neglected. To fix this problem, we invent a hybrid model that 
includes the right fields in moving direction and replaces the transverse field gradients by a 
field strength parameter. In this modification we observe the generation of mono energetic 
electron beams and are able to calculate trapping conditions analytically. 

The next steps of our work will be to go back to the original full gradient model and to 
try to find the envelope approximation functions by other means than the introduction of 
a field strength parameter. If we succeed, trapping cross-sections and trapping conditions 
for the most physical model will be available without thinking about an arbitrarily set field 
strength parameter. 
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